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A strategy to evaluate the distribution of the largest Schmidt eigenvalue 
for entangled random pure states of bipartite systems is proposed. We point 
out that the multiple integral defining the sought quantity for a bipartition 
of sizes N, AI is formally identical (upon simple algebraic manipulations) to 
the one providing the probability density of Landauer conductance in open 
chaotic cavities supporting N and M electronic channels in the two leads. 
Known results about the latter can then be straightforwardly employed in 
the former problem for both systems with broken (/3 = 2) and preserved 
(/3 = 1) time reversal symmetry. The analytical results, yielding a contin- 
uous but not everywhere analytic distribution, arc in excellent agreement 
with numerical simulations. 

PACS numbers: 03.67.Mn, 02.10.Yn, 02.50.Sk, 05.45.Mt, 73.23.-b, 21.10.Ft, 
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1. Introduction 

Consider two sets of n correlated random variables in [0,1], {Aj} and 
{Tj} (z = 1, . . . , n), respectively distributed according to the following joint 
probability densities (jpd): 

(n \ 
Y^X,-l\l[X-ll\Xj-X,\^ (1) 
1=1 / 1=1 j<k 

n 

...,T^) = Ki%l[ Tf n \Tj - Tkf (2) 

i=l j<k 

where Cnf2 and K^^, are known normalization constants and /3 = 1,2. In 
the following two subsection we will provide physical motivations for con- 
sidering such sets, namely set 1 corresponds to the distribution of Schmidt 
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eigenvalues for entangled random pure states in bipartite systems (see sub- 
section 1.1), while set 2 corresponds to the distribution of transmission 
eigenvalues of an open cavity in the chaotic regime (see subsection 1.2). 
Consider now the following statistical quantities: 

• The cumulative distribution Qn{x) = Prob[Aniax < x] of the largest 
member of set 1, Amax = max,{Ai}. By definition, this is given by the 
following n-fold integral: 

Qnix)= [ dAi---(iA„a'i(Ai,...,A„) (3) 
Differentiating Qn{x), one obtciins the pvobability density of A^iax? 

Pn{x) = -J^Qnix) (4) 

• The probability density "Pciy) = Prob[y < G < y + dy] of the quantity 
G = Yll=i which is given by the following n-fold integral ^: 



[ (iri---dT„T2(Ti,...,r„)5 fy- Vr, ) (5) 



Simple algebraic manipulations, summarized in Appendix A lead to the 
following relation between the two quantities above: 



Qn{x) = ^ ) ^ 1/n < X < 1, « = «' 



(6) 



The identity in eq. (6) is the main result of this paper^. Notwithstanding 
its remarkable simplicity, eq. (6) actually permits an exact evaluation of 
Qn{x), the so far unavailable distribution of the largest Schmidt eigenvalue 
for finite n (see subsection 1.1), in terms of 'Pciy) (the probability density 
of Landauer conductance, see subsection 1.2) about which much more is 
known. 

The plan of the paper is as follows. In the next two subsections, we 
give a rather detailed introduction to the physics of entangled random pure 
states in bipartite systems (related to the set 1 above) and about Landauer 
conductance in chaotic mesoscopic cavities supporting a finite number of 



^ One has the foUowing normahzations J^^^ dx Pn{x) = Jg" dy yaiv) = 1- 
^ The bound x >l/n will be discussed in detail later on. 
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electronic channels in the two attached leads (related to the set 2 above). 
In section 2 we exploit the identity (6) to derive analytically the cumula- 
tive distribution and the density of Amax in a few illustrative cases. These 
results are then compared with numerical simulations with excellent agree- 
ment. Eventually we present concluding remarks in section 3 and technical 
developments in the three appendices. 

1.1. Entangled random pure states. 

Entanglement of pure bipartite systems is one of the most active areas 
of research nowadays, due to possible applications to quantum information 
and quantum computation problems [1, 2]. It is also probably the simplest 
setting where well-behaved entanglement quantifiers can be defined, such as 
the von-Neumann or Renyi entropies of either subsystem [2], the so called 
concurrence for two-qubit systems [3] or other entanglement monotones [4, 
5]. 

Typical properties of such states are best addressed by considering ran- 
dom pure states (see e.g. [6] for an excellent review). More precisely, con- 
sider a bipartition of a A^M-dimensional Hilbert space 
:K^P®'Kf\ where we assume without loss of generality that N < M. For 
example, A may be taken as a given system (say a set of spins) living in 
an external environment (e.g., a heat bath) B. A quantum state of the 
composite system can be expanded as a linear combination 

N M 
i=l a=l 

where ji"^) and are two complete basis of and 'K^^^^ respectively. 
The coefficients Xj^a's of this expansion form the entries of a rectangular 
{N X M) matrix X. 

We consider here entangled random pure states This means that: 

1. IV') cannot be expressed as a direct product of two states belonging 
to the two subsystems A and B. 

2. the expansion coefficients random variables drawn from a 
certain probability distribution. 

3. the density matrix of the composite system is simply given by p = 
1-0) ("01 with the constraint Tr[p] = 1, or equivalently (V'lV') = 1- 

More precisely, the density matrix of IV') can then be straightforwardly 
expressed as 

^ = E E 4/3 1 0'^ (8) 
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where the Roman indices i and j run from 1 to iV and the Greek indices a 
and /3 run from 1 to M. 

The reduced density matrix = Tr b[p] is obtained by tracing out the 
environmental degrees of freedom (i.e. those of the subsystem B): 

M 

PA = TrB[p] = Y.{a^\p\a^). (9) 

a=l 

Using the expansion in Eq. (8) one gets 

N M N 

p^=Y. E^^."4«i^^)o-^i = E (10) 

ij=la=l jj=l 

where VFj^'s are the entries of the N x N matrix W = XXK In a similar way, 
one could obtain the reduced density matrix pB = "^^aIp] of the subsystem 
B in terms of the M x M matrix W = X^X. The two matrices W and W 
share the same set of nonzero (positive) real eigenvalues {Ai, A2, . . . , Aat}. 
In the diagonal basis, one can express pA as 

N 

PA = Y.k\Xt){Xt\ (11) 

1=1 

where |A^)'s are the normalized eigenvectors of W = XX^ and similarly for 
Pb- The original composite state |^) in this diagonal basis reads: 

N 
i=l 

Eq. (12) is known as the Schmidt decomposition, and the normalization 
condition (V'l^) = Ij or equivalently Tr[p] = 1, imposes the constraint on 
the sum of eigenvalues, Aj = 1. 

For random pure states, the expansion coefficients in eq. (7) can be 
typically drawn from an unbiased (so called Hilhert- Schmidt) distribution 
(real or complex) 

Prob[X] oc 6 (^T\-(XXt) _ 1^ (13) 

where the Dyson index /3 = 1,2 corresponds respectively to real and com- 
plex X matrices'^. The meaning of eq. (13) is clear: all normalized density 



These two cases in turn correspond to quantum systems whose Hamihonians preserve 
(/3 = 1) or break (/? = 2) time-reversal symmetry. 
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matrices compatible with unitary invariance are sampled with equal proba- 
bility, which corresponds to having minimal a priori information about the 
quantum state under consideration. This in turn induces nontrivial correla- 
tions among the Schmidt eigenvalues (which are now real random variables 
between and 1 whose sum is 1) and makes the investigation of several 
statistical quantities about such states quite interesting. Here we present a 
quick summary of known results: 

• the joint probability density (jpd) of Schmidt eigenvalues, derived by 
Lloyd and Pagels [7], which is precisely given by yi(Ai, . . . , Aat) in eq. 
(1), with n = iV and a = (/3/2)(7V - M + 1) - 1. Note that the delta 
function there guarantees that Tr[/9] = 1. The normalization constant 
in this case reads [8]: 



• the average von Neumann entropy for large iV, M (computed by Page 
[9] for /3 = 2 and extended in [10] to the case /? = 1); 

• the average von Neumann entropy for finite N, M and /3 = 2, conjec- 
tured by Page [9] and independently proven by many researchers soon 
after [11] also in a non-extensive setting [12]; 

• density of Schmidt eigenvalues (one-point function) for finite {N,M), 
derived independently in [13] and [14] for (3 = 2 and in [15] for /3 = 1; 

• universality of eigenvalue correlations for (3 = 2 [IG]; 

• average fidelity between quantum states [17] and distribution of so- 
called G-concurrence [4] for /5 = 2; 

• distribution of so-called purity for small [IS], and phase transitions 
in its Laplace transform for large N [19]; 

• full distribution of Renyi entropies (including large deviation tails), 
computed in [20] for large N = M and all /3s using a Coulomb gas 
method. As a byproduct, the authors also obtain in [20] the average 
and variance of Renyi entropy valid for large = M, and the density 
of Schmidt eigenvalues for all /3's and = cM large; 



N,a={l3/2){N-M+1)-1 



r(MAf/3/2)(r(l + /3/2)) 



n;=o r((M - j)/3/2)r(l + (iV - i)/3/2) 



(14) 



distribution of smallest eigenvalue (related to so-called Demmel con- 
dition number [21]) for /? = 1,2 and finite M = N, derived indepen- 
dently in [22] and [23]. In the latter paper, a conjecture by Znidaric 
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[24] was proven (see also [2-"')] for an extension of these results to the 
case iV / M); 

• distribution of largest eigenvalue for large N = M and all /5s [20], 
including small and large deviation laws. Typical fluctuations around 
the mean ~ A/N (once properly scaled) are found to follow the Tracy- 
Widom distribution (see also [26] for a related result). No results seem 
to be available for the case of finite N, M. 

Given the current interest in the distribution of extreme Schmidt eigenval- 
ues, the reader may on the one hand wonder whether they really encode 
useful information, and on the other why the largest eigenvalue distribution 
for finite A^, M is much harder to obtain via the same strategy used for the 
smallest one [23]. 

In order to answer the first question, first note that due to the constraint 
YliLi = ^ aiid the fact that all eigenvalues are nonnegative, it follows 
that^ < Ajnax ^ 1 and < Amin ^ Now consider the following 

limiting situations. Suppose that the largest eigenvalue Amax = niaxjjAi} 
takes its maximum allowed value 1. Then it follows immediately that all 
the remaining {N — 1) eigenvalues must be identically 0. In this situation 
eq. (12) tells us that {tp) is fully unentangled (completely separable). On 
the other hand, if Amax = 1/-^ (i-e., it takes its lowest allowed value), all 
the eigenvalues must have the same value, Aj = 1/A^ for all i. In this 
case, the pure state {ip) is maximally entangled, as this state maximizes 
the von Neumann entropy S'vn = — ^^^i = In(A^). In other words, 

the knowledge of the largest eigenvalue distribution really provides essential 
information about how entangled a random pure state is. 

A discussion about the asymmetry in the treatment of smallest and 
largest Schmidt eigenvalues is included in Appendix B. 

1.2. Landauer conductance in open cavities 

Consider a cavity of submicron dimensions etched in a semiconductor 
and connected to the external world by two leads supporting M and N 
electronic channels. It is well established that the electrical current flowing 
through such a cavity when brought out of equilibrium by an applied exter- 
nal voltage presents time-dependent fluctuations which persist down to zero 
temperature [27] and are thus associated with the granularity of the electron 

" For /3 = 2, (Amiii> = exactly for all N = M, while for /? = 1 (A,„in) ~ Ci/Af^ for 

large N = M, where the constant ci is precisely known [2.?]. 

^ This means that both the smallest and the largest eigenvalue distributions have com- 
pact supports and justifies the bound x > 1/n in eq. (6). 
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charge e. Typical features observed in experiments include weak localiza- 
tion [28], universality in conductance fluctuations [29] and constant Fano 
factor [30]. The Landauer-Biittiker scattering approach [27, 31, 32] is rather 
successful in describing the statistics of quantum transport: it amounts to 
relating the wave function coefficients of the incoming and outgoing elec- 
trons through the unitary scattering matrix S (2A''o x 2A''o, if A'^o = N + M): 

where the transmission (t,t') and reflection (r, r') blocks are submatrices 
encoding the transmission and reflection coefficients among different chan- 
nels^. Experimental observables can be extracted from the eigenvalues of 
the hermitian transport matrix T = tt^: for example, the dimensionless 
conductance and the shot noise are given respectively by G = Tr(T) [31] 
and P = Tr[T(l - T)] [:]:■]. 

Random Matrix Theory (RMT) has been very successful in describing 
the statistics of universal fluctuations in such systems, when the correspond- 
ing classical dynamics is chaotic: the scattering matrix S is drawn from a 
suitable ensemble of random matrices, with the overall constraint of unitar- 
ity [34, 35, 36]. A maximum entropy approach (under the assumption of 
ballistic point contacts [27]) forces the probability distribution of S to be 
uniform within the unitary group, i.e. S belongs to one of Dyson's Circular 
Ensembles [37, 38]. 

Prom the uniformity of the probability density of § within the unitary 
group, the jpd of the transmission eigenvalues {Tj} of the matrix T, from 
which the statistics of interesting experimental quantities could be in prin- 
ciple derived, is precisely given by eq. (2), with n = min(iV, M) and 
a' = f (|iV - M\ + 1) - 1 [27, 36, 39]. There, the Dyson index /? char- 
acterizes different symmetry classes (/3 = 1,2 according to the presence or 
absence of time-reversal symmetry). The eigenvalues Tj are thus correlated 
real random variables between and 1. The normalization constant K^^^, 
is explicitly known from the celebrated Selberg's integral as: 

_ TT r(l + /3/2)r(2 + a' + {/3/2){n + j - 1)) 
n,a' lir(l + (/3/2)(l+i))r(l + a' + /3j72)r(l + /3j72) ^ > 

From (2), in principle the statistics of many observables of interest can 
be calculated. In particular, we focus here on the dimensionless Landauer 
conductance G = Tr(T) = Y17=i which satisfies the bounds < G < n. 



(t, t') are respectively of size N x M and Af x A'^, while (r, r') are of size M x M and 
N X N. 
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Mean and variance of G and other quantities have been known for a long 
time when N,M are large [27, 40, 41], and recently also for finite N,M 
[41, 42, 43, 44]. 

Stimulated by some recent experimental progresses [45], which made 
eventually possible to explore the full distribution of the conductance (and 
not just its mean and variance), a lot of effort has been devoted to its 
theoretical characterization. An explicit expression was first obtained for 
n = 1,2 [46, 47, 48], while more results were available in the case of quasi 
one-dimensional wires [49] and 3D insulators [50]. For the shot noise, the 
full distribution was known only for N = M = 1 [51]. Very recently, Som- 
mers et al. [52] announced two formulas for the distribution of conductance 
and shot noise, valid at arbitrary number of open channels and for any /3, 
which are based on Fourier expansions. In [53], a complete and systematic 
approach to the full statistics of conductance and shot noise was brought 
forward. This is based on symmetric function expansions, and is valid for 
(5 = 1,2 and arbitrary M,N. In addition, the authors of [53] were able 
to provide general formulae in terms of determinants or Pfafhans for the 
full probability distribution of conductance and shot noise at quantized /3 
and general N,M (see below). In [54], the integrable theory of quantum 
transport in chaotic cavities (based on Painleve transcendent) for (3 = 2 was 
formulated, and recursion formulae for the efficient computation of conduc- 
tance and shot noise cumulants have been derived. In two recent publica- 
tions [55], the distributions of conductance, shot noise and integer moments 
were computed in the limit = M ^ 1 using a Coulomb gas method, and 
long power-law tails were detected in the distributions^, a result confirmed 
by extensive numerical simulations. In [56], exact results (basically equiv- 
alent to the general formulae in [53]) for the Laplace transform 7g{s) of 
the conductance distribution for any M and all /3s are formally given as 
follows: 



/3 = 1: 



7g{s) oc Fi[^^ll{s)]j^k=o,...,N^i for N even, and Tg{s) oc 



Pf 







for N odd, where F{[A] is the 



j,k=0,...,N-l 



Pfaffian of the even-dimensional antisymmetric matrix A [37] and the 
proportionality constants are known explicitly. The arrays ^^^1 (s) and 



^ A careful asymptotic analysis of formulae in [54] leads to the same findings. In par- 
ticular, there are no long exponential tails in the distributions, as originally claimed 
in [54]. 
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{s) are given by: 



= C r sgn(x - 2/)e-^^e-^^x"'+^y°'+'=dx dy (17) 
Jo Jo 

r-l 

^f{s)= / e-'^x'^'+^dx (18) 
Jo 

• (3 = 2: in this case, we have a representation in terms of a Hankel 
determinant 

7g{s) oc det[^'g(s)]i,fc=o,...,iV-i (19) 

where: 

^fl{s)= / e-^^x^'+^+'^dx (20) 

JO 

that was first derived in [54]. 

Exphcit inversions of the Laplace transforms above are always possible on 
a case-by-case basis, and a catalogue of such evaluations for few interesting 
cases is provided in ref. [56]. In next section and in Appendix C, we will 
combine such explicit formulas with the identity in eq. (6) to illustrate the 
validity of our approach to the distribution of largest Schmidt eigenvalue. 

2. Two applications of the main identity 

For illustrative purposes, we consider two applications of the main iden- 
tity eq. (6). A third one is discussed in great detail in appendix C. 

1. iV = 3, M = 4, /3 = 1: in this case, the probability density of Landauer 
conductance has been derived explicitly in [56] as: 

'^Giy) = l[y'-{y- ifiv^ - 122/ + 51)0(y - 1) - (y - 2)^ 

x(y2 + 6y + 24)0(y-2)]0(3-y) (21) 
where 9{x) is Heaviside step function. 

2. N = 3, M = 5, (3 = 1: again, the probability density of Landauer con- 
ductance in this case has been derived explicitly in [56] as: 

' ^y'^/' for < y < 1 

2^ [3003y5 - 21021/ + 55770^^ 
"^ciy) = I -70070y2 + 42315y - 9933 

-32(y - 2)^^/2 (2y3 + 14^2 + 63y + 231) 
xe{y-2)] forl<y<3 

(22) 
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Note that: 

• since the conductance density has always a compact support [0, n] 
(where n = min(A^, M)), it follows immediately from eq. (6) that the 
cumulative distribution of the largest Schmidt eigenvalue Qn{x) and 
its density Pn{x) have compact support 1/??, < x < 1, as expected. 

• since the conductance density is known to be continuous but not ev- 
erywhere analytic [52, 53] (i.e. it displays 'critical' points at which 
higher derivatives are discontinuous), the cumulative distribution and 
the density of the largest Schmidt eigenvalue enjoy this property too. 

In fig. 1 we plot the cumulative distributions Qn{x) corresponding to 
the two cases above. The curves are obtained via eq. (6), where ydv) 
is respectively given by (21) and (22). The distributions are increasing 
functions of the argument x, as they should, and such that (5n(l) = 1- 
Differentiation of the analytical formulas provide the density of the largest 
eigenvalue, Pn{x), which is plotted in figs. 2 and 3 along with numerical 
simulations. These are obtained as follows [8, 57]: 

1. we generate k ~ 10^, 10^ real Gaussian M x N matrices X. 

2. for each instance we construct the Wishart matrix W = X'^X. 

3. we diagonalize W and collect its real and non- negative eigenvalues 
{Ai, . . . , Atv}. 

4. we define a new set of variables < Aj < 1 as Aj = Aj/^^-^ Aj, for 
i = 1, . . . ,N. The set of variables Aj is guaranteed to be sampled 
according to the measure (1). 

5. we construct a normalized histogram of Amax = niaxj{Aj}. 

The agreement between theory and simulations is excellent. 

We also evaluated exactly the average of (Amax) = fi/n^^ ^ Pn{x) for 
the two cases above and found: 

25 

(Arnax)l = 1^ - 0.694444... (23) 
3d 

(Amax)2 = ^(378 + 89^/3) « 0.656978... (24) 

oiU 

However, a general formula for (Amax) is still elusive (unlike (Amin) [23]) and 
may well attract further researches (see also Appendix C). 
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e„(x) 




0.2 0.4 0.6 0.8 1.0 

Fig. 1. Cumulative distribution Qn{x) of the largest Schmidt eigenvalue via eq. (6) 
for 71 = = 3, 7\/ = 4, /3 = 1 (violet) and n = = 3, il/ = 5, /3 = 1 (blue). This 
distribution is identically zero for x <l/n and Q„(l) = 1 as it should. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8, 0.9 1 



Fig. 2. Density of the largest Schmidt eigenvalue p„(x) via eqs. (4) and (6) for 
n = N = 3, M = 4, [3 = 1 (solid black curve), along with numerical diagonalization 
of K = 8 • 10* samples (red dots), see main text for the algorithm. 



3. Conclusions 

We have presented an exact identity relating two statistical quantities 
which arise in different contexts: the cumulative distribution of the largest 
Schmidt eigenvalue for entangled random pure states in bipartite systems 
of sizes (M, N) and the probability density of Landauer conductance in 
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Fig. 3. Density of the largest Schmidt eigenvalue Pn{x) via eqs. (4) and (6) for 
n = N = 3,M = 5,f3 = l (solid black curve), along with numerical diagonalization 
of K = 10^ samples (red dots), see main text for the algorithm. 



chaotic cavities supporting N and M electronic channels in the two exter- 
nal leads. Recent analytical results for the latter are exploited to derive 
(so far unavailable) exact formulas for the former quantity at finite N, M 
(while large N, M results are already available [20] ) , which is of interest 
in order to quantify the degree of entanglement of random pure states. A 
detailed introduction to the physics involved has been provided, along with 
a precise discussion of the asymmetry in the treatment of the smallest and 
largest Schmidt eigenvalue distributions. A general formula for (Amax 

), the 

average of the largest eigenvalue, valid for arbitrary N, M is unfortunately 
still lacking and certainly deserves further investigations. 
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manuscript and many helpful advices. 



Appendix A 

Derivation of main identity. 
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Consider the cumulative distribution Qn{x) of Amax = niaxj{Aj} (eq. 
(3)): 

// n \ " 

dAi • • • d\n5 [Y.^^-A\{K\{\>^j- Afel^ (A.l) 
^J'^l" \j=l / j=l j<k 

A change of variables Aj = xTi leads to: 

x/ dTi---dTn5[xy^Ti-l]\[T^T[\Tj-n\^ (A.2) 

■^[0.1]" V ^ /r=i 7<fc 

Using the property of delta functions 5{'^Y) = (5(y)/|7|, we obtain straigh- 
forwardly eq. (6). 

Appendix B 

Asymmetry in the treatment of smallest and largest Schmidt eigenvalues 

Consider the cumulative distribution Qn{x) of the smallest member of 
set 1, Ainin = niinj{Aj} (i.e. the distribution of the smallest Schmidt eigen- 
value for entangled random pure states). By definition it is given by: 

0n(x) = / dAi • • • d\n5 A, - 1 n A^ n l^i - ^kf (B.l) 



i=l / i=l j<k 



where the upper limit of integration can be safely extended up to oo in view 
of the unit norm constraint. In the case a = and /3 = 2, which we focus 
on here for illustrative purposes, the evaluation of this multiple integral 
proceeds via the auxiliary function [2,3]: 



f J / dAi • • • d\n5 ( E - t I n \\j - Afc|2 (B.2) 



such that Qnix) = G.„(3;,l). Next, one takes the Laplace transform of 
@n{x,t): 



f oo 

dtQr 



'0 



{x,t)e-'' = Ci'i [ dXi ■ • -dA^e-^^-i^^ JJ \Xj - \k\^ (B.3) 
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and in the r.h.s performs a linear shift yj = s{Xi — x) (this is the crucial 
technical step), to get: 



/■oo —slMx r 

/ dt@nix,t)e-'' = -^Cl^l / dy^--- dy^.e-^'^^y^ TT \y, - y^l^ 

Jo S'^ J[0^^]n J-^-^ 

' V ' 

ocZwL 

(B.4) 

Note that, thanks to the linear shift, the dependence on the Laplace variable 
s has been entirely transferred outside the n-fold integral: this is now pro- 
portional to the partition function 2.wL of an associated Wishart-Laguerre 
(WL) ensemble of random covariance matrices [37, 58] of the form W = X^X, 
where X is a Gaussian rectangular matrix with real or complex entries. The 
joint distribution of the n nonnegative eigenvalues of W is known [59] 

T(WL)(^^, . . . , = y^m e-l EIL. 'flXtll |A, - X,f (B.5) 

i=l j<k 

where JSf^^a is a known normalization constant. Therefore, the jpd of Schmidt 
eigenvalues (1) can be seen as a fixed-trace (microcanonical) version of the 
Wishart-Laguerre (canonical) ensemble^. 

Now, in the case of the cumulative distribution of the largest eigenvalue 
Qn{x), the integrals on the r.h.s. run over [0,x] instead of [x,oo], making 
the aforementioned linear shift less useful. One could keep pursuing the 
Laplace-transform route (introducing an auxiliary function Qn{x,t)) with 
the change of variables sXi = yi obtaining: 

/ dtQr,{x,t)e-'' = ^ dyi---fiy„e-^-i2'^rr|y,-yfc|MB.6) 



o^-BwL[(sa;,oo)] 



but the integral on the r.h.s. does not permit this time a friendly Laplace- 
inversion (see however next appendix). Indeed, this integral is readily recog- 
nized as proportional to £'wl[(s2;, cxd)], where Swl[(o, 6)] is the gap proba- 
bility for a Wishart-Laguerre ensemble, i.e. the probability that the interval 
(a, b) on the real axis is free of eigenvalues. This quantity is exactly known 
in terms of Painleve V [62], with the consequence that an explicit Laplace 
inversion formula is not available to date. 



Note that the presence of a fixed-trace constraint has crucial consequences on the 
spectral properties of random matrix ensembles [60, 61]. 
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In summary, the asymmetry in the treatment of the smahest and largest 
Schmidt eigenvalues arises because the linear shift that works in the for- 
mer case fails in the latter, and this fact calls for the alternative approach 
developed in this paper. 



Appendix C 

A third application of main identity for (3 = 2 and N = M 

In this appendix we discuss in more detail a third application of the main 
identity eq. (6). In the case (3 = 2 and M = N, the Hankel determinant 
representation in eq. (19) actually allows a more systematic (and more 
easily automatized) treatment of the cumulative distribution of the largest 
Schmidt eigenvalue. Following the Laplace transform route outlined in the 
previous appendix, one can easily write down the following equation: 



r dtQ.^{x, t)e-''/'- = Ci^l x"' / dyi--- dy^e-' ^"-1 ^» TT \y, - y^f 



n! dct[(-a,y+'=:fi(s)]j.fe=o,....n-i 

(C.l) 

where we have used the Hankel determinant representation for the n-fold 
integral on the r.h.s. provided in [54] (fully equivalent to (19)), with 3"i(s) = 
{l-e-^/s. 

A change of variable t = xr on the l.h.s. followed by a (formal) Laplace 
inversion, leads to the final formula for Qn{x) in this case^: 

Qn{x) = Cj^l n\ x"'^i£-i [det[(-a,)^'+^'Ji(s)],-fc=o,...,n-i] (V^) (C.2) 

where L~^[f{s)]{t) is the inverse Laplace transform of f{s) with parameter 
t. The r.h.s. of (C.2) can be systematically evaluated in Mathematica® 
and its derivative {pn{x)) has been plotted in fig. 4 for n = N = M = 3,4. 
Also the average of the largest eigenvalue (Amax) = Ji/n'^^ ^ Pn{x) can 

be evaluated exactly in Mathematica® for a given n, and we provide a 
few evaluations in the following table. Unlike the average of the smallest 
eigenvalue, which in the same circumstances (/3 = 2 and M = N) has the 
attractively simple form (Amin) = 1/-/V^ exactly for all N [23], the situation 
for the largest eigenvalue seems more complicated and it is much harder 
to even conjecture a possible formula for the sought average valid for all 



^ One could have started directly from eq. (6) and would have been 
led to the very same eq. (C2) upon noticing that Vci^/x) = 

{n\/K'^l)L-^ [det[(-9.)^+'--Ji(.)],,fe.o,...,n-i] {1/x) 
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eq. (C.2) for n = = A/ = 3 (blue) and n = N = M = i (violet) and /? = 2. 



= M. We leave this as a challenging open problem, noticing en passant 
that the asymptotic value (Amax) ~ derived in [2(1] is approached very 
slowly as N increases [15]. This fact makes the investigation oi finite (small) 
N results very much called for. 



N = M 


(Amax) (exact) 


(Amax) (approx.) 


2 


7 

8 


0.875 


3 


313 
432 


0.724537 


4 


1367807 
2239488 


0.610768 


5 


4581882694877 
8707129344000 


0.526222 


6 


225128892964655720357665283 
487487792008396800000000000 


0.461814 
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